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ABSTRACT 


Ultraviolet  absorption  was  used  to  measure  the  rate  of  formation  of  difluoro- 
methylene  (CF2)  from  decomposed  CF31,  C2F6,  and  CF4  in  excess  argon  behind  shock 
waves.  In  some  experiments  pure  CF4  was  shocked.  Data  were  taken  over  a  temp¬ 
erature  range  from  1700°  to  3000°  K  at  total  concentrations  between  2(10”5)  and 
5(10“^)  mole-cm”3.  A  chemical,  nonequilibiium  shock-tube  computer  program  was 
developed  to  analyze  the  CF2  kinetic  profiles.  By  curve  fitting  the  data,  rate 
constants  for  a  number  of  fluorocarbon  reactions  are  obtained. 

Unclassified  report 

A  copy  of  this  report  is  available  at  Reports  Distribution  Center,  Room  1126 
and  will  be  held  for  a  period  of  three  weeks. 
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INTRODUCTION 


The  role  of  fluorine  in  the  high  temperature  oxidation  chemistry  of  tetrafluoro- 
ethylene  (C2F4)  is  little  understood.  For  example,  among  the  reaction  products 
of  C2F4  and  its  polymer  with  oxygen,  there  have  been  identified  CF4  and  CF3 
species. (1*2)  it  was  suspected  that  these  molecules  resulted  from  the  reaction 
between  fluorine  and  the  dif luoromethylene  radical  (CF2).  Support  of  this 
hypothesis  was  noted  when  shock  heated  mixtures  of  CF3I,  and  CF4  in  argon 

diluent  as  well  as  pure  CF4  were  found  to  decompose  to  CF2.  Experiments  were 
carried  out  to  detect  the  CF2  radical  by  its  absorption  in  the  ultraviolet.  In 
the  present  paper,  measured  kinetic  profiles  of  this  species  are  analyzed  with  a 
computer  program.  The  computer  program  solves  simultaneously  the  Hugoniot  rela¬ 
tions  (mass,  momentum,  and  energy)  and  the  chemical  rate  equations  for  a  chemically 
unrelaxed  shock  wave.  The  calculated  rate  profiles  of  CF2  are  compared  with  ex¬ 
periment.  Some  of  the  rate  constants  employed  in  the  analysis  were  taken  from 
literature  sources,  some  were  determined  empirically  and  others  were  estimated 
from  theory  or  deduced  from  JANAF  equilibrium  constants. 
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EXPERIMENTAL  MEASUREMENTS 


The  shock-tube  apparatus  has  been  described  previously. It  consists  of  1.5- 
inch  inside  diameter  stainless  steel  sections  provided  with  lij’.ht  screens  for 
shock  velocity  measurements  and  an  optical  station  for  ultraviolet  absorption 
spectroscopy.  Gas  mixtures  of  CF4  (Matheson,  95  percent  purity),  CF3I  (Air  Products, 
98  percent  purity),  CoFg  (Air  Products,  99.9  percent)  in  excess  argon  (Matheson, 
99.999  percent  purity)  were  prepared  for  study.  A  few  tests  were  conducted  in 
undiluted  CF4.  The  CF2  radicals  in  the  shock-heated  reaction  matures  were  fol¬ 
lowed  photometrically  in  absorption  at  266oX.  The  absorption  coefficient, 

1.15  x  106  (±10  percent)  cubic  centimeter  per  mole-centimeter,  was  used  to  calcu¬ 
late  concentration  from  Beer’s  law. w)  Oscillogram  records  shown  in  Figure  1  are 
typical  of  the  CF2  absorption  from  these  gas  mixtures.  The  CF3I  and  C2F5  molecules 
were  studied  since  CF3  radicals  are  immediately  formed  behind  the  shock  wave, 
making  the  rate  determining  step  for  CF2  appearance  the  dissociation  of  CF3. 

With  CF4,  the  complete  fluorocarbon  chemistry  involving  CF4,  CF3  and  CF2  is  given. 

By  studying  the  decomposition  rates  leading  to  CF2,  the  kinetics  of  the  fluorine 
oxidation  reactions  may  thus  be  ascertained  from  detail  balancing  with  appropriate 
equilibrium  constants. 
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NONEQUILIBRIUM  SHOCK  TUBE  PROGRAM 


The  equations  governing  a  moving  one-dimensional  inviscid  fluid  in  a  shock  tube 
are  the  so-called  Rankine-Hugoniot(^)  relations  expressing  conservation  of  mass, 
momentum,  and  energy. 


d(pu)/dt 

=  0 

(mass) 

(1) 

pu  d ( u + 

—  )/dt  =  0 

(momentum) 

(2) 

pu 

pu  d(h  +  1/2  )/dt  m  0 

(energy) 

(3) 

where  p  is  the  density,  u  is  the  flow  velocity,  p  is  pressure  and  h  is  the  en¬ 
thalpy  per  unit  mass.  For  a  multi-component  gas  h  is  defined  by 


2 


(4) 


The  quantity  c(-  is  the  concentration  of  the  ic^  species  in  mole  per  cubic  centi¬ 
meter,  hi  is  the  molar  enthaipy  and  Mi  the  gram  molecular  weight.  The  ideal  gas 
equation 


P 


P 


RT 


(5) 


is  used  to  relate  the  thermodynamic  temperature  T  to  pressure  and  density  where 
R  is  the  universal  gas  constant  per  mole.  The  molar  enthalpy  is  then  expressed 
in  terms  of  temperature  by  a  power  series, 


h-  - 


-  h  T~‘ + 


(6) 


A  curve  fit  of  the  molar  enthalpy  of  each  species  was  generated  from  JANAF  thermo¬ 
chemical  data. (5)  The  variation  in  species  concentration  resulting  from  chemical 
reaction  and  volumetric  change  behind  the  shock  wave  is  given  by 


dci  ci  dp 

dt  p  dt 


1  nvkj  ek  Aq  T  i  *xp(Efj/RT) 
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where  the  positive  term  on  the  right  of  the  equality  reflects  the  sum  of  the 
chemical  rates  leading  to  species  production  and  the  negative  term,  the  rates  for 
removal.  Subscripts  /  and  r  denote  the  forward  and  reverse  processes,  ui;  and 
vq  are  the  stoichiometric  coefficients  for  the  and  ith  species  in  the 
reaction;  Aj ,  nj  and  Ej  are  the  pre-exponential  factor,  temperature  exponent,  and 
activation  energy  of  the  rate  constant  for  that  reaction.  To  solve  this  set  of 
differential  equations  for  a  multicomponent,  chemically  relaxing  shock  wave,  a 
Runge-Kutta  integration  scheme^)  and  an  IBM  360  computer  were  employed.  The  in¬ 
puts  to  the  calculation  are  the  shock  velocity,  state  of  the  gas  ahead  of  the 
shock  wave,  the  chemical  reactions  to  be  considered,  the  forward  and  reverse  rate 
constants,  and  coefficients  of  the  enthalpy  curve-fits.  Typical  machine  times 
for  solving  the  fluorocarbon  chemical  kinetics  problems  *an  between  1  and  2 
minutes.  Before  using  the  nonequilibrium  shock- tube  program  to  analyze  the  present 
data,  a  progrc-1  check  was  made  against  an  air  chemistry,  stream  tube  program 
developed  by  !'  namin  and  O'Brien.^)  Figure  2  shows  that  calculations  from  both 
programs  for  electron  density  behind  a  shock  wave  into  air  agree  numerically  to 
within  3  percent. 


CHEMICAL  MECHANISMS  AND  KINETICS 


The  reactions  considered  in  this  study  are  shown  in  Table  1. 

CF3I  and  ^2^6  Decomposition:  A  bimolecular  rite  constant  for  R1  was  estimated 

from  classical  collision  theory assuming  energy  contributions  from  half  the 
number  of  vibrational  modes  and  a  steric  factor  of  0.1.  The  dissociation  rate 
constant  of  R2  was  obtained  from  the  JANAF  equilibrium  constant  and  the  rate 
constant  of  Ayscough^'  for  trifluoromethyl  radical  recombination.  Calculations 
based  on  these  rate  constants  (Figures  3  and  4)  showed  that  within  a  few  micro¬ 
second?  after  shock  passage  CF3I  was  dissociated  by  several  orders  of  magnitude 
and  that  was  about  80  percent  dissociated  with  R2  and  R2'  in  equilibrium. 

Under  these  conditions,  the  apparent  rate  constant  for  CF2  formation  by  CF3 
dissociation  (R3)  is  given  by 


d[CF2) 

- - -  «k3(expt’t)  \[CF3]0-(I+y)[CF2)\[Ar)  (8) 

P21  dtf 


where  y  is  zero  if  R4'  «  R3  and  one  if  F4'»RJt  i.e.,  for  each  CF2  produced  one 
CF3  radical  is  removed  when  fluorine  attachment  is  slow  and  two  CF3  radicals  are 
removed  when  fluorine  attachment  is  fast.  The  concentration  of  CF3  initially 
behind  the  shock  wave  is  [CF3  The  density  ratio  P21  is  multiplied  in 

Equation  (8)  to  reference  the  event  to  laboratory  time  (*£ ) .  In  calculating  k3 
(expt'l)  by  means  of  Equation  (8),  y  wa3  taken  to  be  unity  for  the  lower  tempera¬ 
ture  data,  and  was  assumed  zero  for  the  higher  temperature  data.  Calculations 
(Figures  3  and  4)  indicated  that  when  analyzing  the  low  temperature  data  after 
50  percent  decomposition,  the  above  approximation  was  good  to  within  90  percent, 
.'mother  feature  brought  out  by  the  calculations  was  that  after  CF^I  ar.d  C2F^ 
chemical  relaxation,  the  shock  temperature  profile  was  nearly  constant,  so  that 
little  temperature  uncertainty  was  introduced  in  evaluating  the  rate  constant 
for  R3  during  reaction.  With  the  higher  temperature  data,  k3  (expt'l)  was 
determined  from  the  initial  CF2  absorption  slope.  As  seen  in  Figure  5,  the 
rather  good  agreement  obtained  within  the  scatter  limits  of  k3  (expt'l)  in  the 
CF3X  and  l^F^  experiments  tends  to  uphold  the  data  reduction  method.  Examples 
of  the  CF3I  and  C2Fg  study  are  given  in  Table  2. 

CF^  Decomposition:  In  the  temperature  range  2200°  to  3000°  K,  the  CF3  dissocia¬ 
tion  rate  is  three  orders  of  magnitude  faster  than  the  rate  of  appearance  of  CF2 
in  the  oscillogram  records  of  CF^  decomposition.  Hence  these  observations  sug¬ 
gested  that  the  rate  determining  step  for  CF2  formation  starting  with  CF^  was 
R4.  Calculations  (Figure  6)  showed  that  under  the  experimental  conditions  chosen 
the  CF3  co^nt ration  would  be  much  less  than  that  of  CF2  and  CF4  in  the  reaction 
mixture,  allowing  an  apparent  rate  constant  for  R4  to  be  evaluated  from  the 
initial  CF2  slope  measurements  according  to 

d[CF2] 

—7—  »  k4  (expt'l)  l[CF4]0  ~  [CF2]\  [ Ar ].  (9) 

P21  dtl 
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Examples  of  the  experiments  are  shown  in  Table  3.  The  temperature  dependence  of 
the  apparent  rate  constants  obtained  from  1:100  and  2:100  CF^Argon  gas  mixtures 
is  shown  in  Figure  7.  The  kinetic  profile  of  CF2  measured  in  a  pure  CF^ 
run  (Figure  8)  is  compared  to  that  calculated  from  the  nonequiliblrum  shocktube 
program  using  the  rate  constants  in  Table  1.  .  The  results  indicate  that  the 
collision  efficiencies  of  CF4  and  argon  appear  to  be  quite  similar.  Another 
aspect  of  the  pure  CF4  experiment  is  that  the  intermediate  CF  chemistry  becomes 
important  and  serves  as  a  check  on  its  rate  constants  in  the  calculations.  The 
rate  constants  for  R5  and  R5'  were  obtained  in  an  earlier  study  of  the  thermal 
decomposition  of  the  difluoromethylene  radical^10'. 

Finally,  for  a  number  of  shock  conditions,  measured  and  calculated  kinetic  pro¬ 
files  cf  CF2  from  decomposition  of  the  three  fluorocarbon  molecules  investigated 
are  shown  in  Figure  9.  It  is  seen  that  the  rate  constants  given  in  Table  1 
would  appear  to  be  adequate  in  describing  the  experimental  data  over  the  tempera¬ 
ture  range  of  the  study. 
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DISCUSSION 


The  heat-of-r°action  for  R3  was  obtained  from  JANAF  thermochemical  data  and  used 
to  fit  the  experimental  rate  constants  of  the.CF^I  and  C2*(>  data.  The  large 
temperature  exponent  (t“9.04)  determined  appears  to  be  Inconsistent  with  the 
classical  collision  treatment  for  a  polyatomic  like  CF^.  The  maximum  number  of 
internal  degrees  of  freedom  not  counting  rotations  which  could  contribute  energy 
for  reaction  would  be  6.  Rotational  and  translational  participation,  if  included, 
would  increase  the  temperature  exponent  to  only  (T”8.5).  jt  was  found  that  in 
calculating  the  CF2  rate  profiles,  good  agreement  with  the  measurements  could  be 
obtained  with  the  JANAF  equilibrium  constant  multiplied  by  16.  With  the  equili¬ 
brium  pr^-exponential  factor  kept  constant,  a  factor  of  16  gives  a  heat-of-reaction 
of  81.2  kcal/mole  at  2000°  K  (average  temperature  range  of  data).  By  refitting 
the  rate  constants  of  R3  with  the  revised  heat-of-reaction,  a  temperature  ex¬ 
ponent  of  (T-6.21)  is  obtained,  which  is  more  in  accord  with  theory. 

Besides  the  reactions  shown  in  Table  1,  the  rate  data  were  also  analyzed  by  add¬ 
ing  the  reactions 


R6,6’ 

CF4  +  F  * 

cf3  +  f2 

R7,7f 

CF3  +  F  £ 

cf2  +  F2 

R8,8' 

F2  +  M  ^ 

F  +  F  +  M 

Calculations  based  on  rate  constants  from  classical  collision  theory  and  JANAF 
equilibrium  constants  showed  that  R6,6'  and  R7,7*  were  about  three  orders  of 
magnitude  slower  than  the  decomposition  reactions  and  had  negligible  effect  on 
the  measured  CF2  kinetic  profiles. 
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APPENDIX  A 


Formulation  of  the  Nonequilibrium  Shock-Tube  Program 

The  nonequilibrium  shock  tube  program  solves  simultaneously  the  Rankine-riugonio t 
conservation  relations  and  the  chemical  rate  equations  expressing  variation  In 
s  ides  concentration.  The  program  is  organized  to  Integrate  the  rate  equations 
stepwise,  an  iteration  scheme  to  satisfy  the  conservation  relations  being  repeated 
at  each  step.  It  differs  from  the  average  ordinary  differential  equation  program 
chiefly  in  that  the  equations  themselves  are  not  fixed.  Instead,  only  rules  for 
generating  the  equations  are  fixed. 

The  input  quantities  include  the  shock  velocity;  the  temperature,  pressure,  and 
species  mole  fractions  ahead  of  the  shock;  the  chemical  reactions  to  be  considered; 
the  preexponential  factor,  temperature  exponent,  and  activation  energy  of  the  rate 
constant  for  each  reaction;  and  coefficients  of  the  enthalpy  curve  fits.  A  library 
tape  of  species  data  (molecular  weights  and  enthalpy  fit  coefficients)  may  be 
searched  for  quantities  omitted  from  the  regular  program  input.  The  chemical 
reactions  are  entered  directly  in  equation  form,  a  maximum  of  four  species  allowed 
on  each  side  of  an  equation,  with  each  species  allotted  a  maximum  of  eight  alpha¬ 
numeric  characters.  For  example,  consider  the  simple  system 

C2F4  +  M  *  CF2  +  CF2  +  M 

CF2  +  CF2  +  M  -  C2F4  +  M 

CF2  +  M  «  CF  +  F  +  M 

CF  +  F  +  M  -  CF2  +  M 

where  M  will  be  taken  by  convention  to  refer  to  the  sum  of  all  species  concentra¬ 
tions  . 

The  set  of  equations  is  analyzed  by  the  program,  the  species  in  each  equation 
matched  against  an  input  list  of  all  species  occurring  in  the  problem,  and  a 
matrix  M  of  indices  set  up  to  represent  the  original  equations.  If,  in  the  ex¬ 
ample  above,  the  species  are  arbitrarily  ordered, 

C2F4  CF2  CF  F 

1  2  3  4 


and  M  is  assigned  the  index  5, 


M  - 


15  0  0 
2  2  5  0 

2  5  0  0 

3  4  5  0 


2  2 

1  5 
3  4 

2  5 


the  4x8  matrix  generated  by  the  program  is 
5  0 
0  0 
5  0 
0  0 
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The  rules  for  generating  the  set  of  differential  equations  are  then  as  follows: 

For  the  * **  species  the  change  in  the  concentration  C ,•  (t)  as  a  function  of  time 
is  defined  as 

"R  4 

i  ciM  ■  2  v  n  \K  *- '•■■■•  n* 

j = 1  1 


where 

Ns  is  the  number  of  species 
Nr  is  the  number  of  reactions 

Rj  is  the  rate  constant  corresponding  to  the  equation 

-n  if  the  species  occurs  n  times  on  the  left 
side  of  equation  j 

if  the  itfl  species  does  not  occur  on  the  left 
o  side  of  equation  j 

+  n  if  the  itfl  species  occurs  n  times  on  the  right 

side  of  equation  j 

if  the  species  does  not  occur  on  the  right 
o  side  of  equation  j 

and  by  definition 

c„  -  1 


Again,  for  the  example,  the  equations  are 


4 

( -  Rj  Cj  + 

i=l 
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4 

( 2  R,  Cj  -  2  R2  C22  -  R3  C2  +  R4  C3  C4)  Ci 

i  *  / 

4 

(R3C2  -  r4  c3  C4)  2  ci 
i -l 

4 

—  c4  =  (R3  c2  -  r4  c3  C4 )  ^  ^  C; 

*  - 1 


To  find  initial  conditions  for  the  system  of  differential  equations,  the  con¬ 
servation  relations  are  first  solved  for  values  of  the  temperature  T2  (0)  and 
density  P2(0)  .  A  false  position  iteration  for  temperature  is  used  in  the  solution. 
Then,  since  the  mole  fractions  Xi  (t)  are  known  (input)  quantities  at  time  0  , 
the  concentrations  may  be  found  from  the  relation 

Xt(t)  p2(t) 

X;  (t) 

»=  *  .L 

where  is  the  molecular  weight  of  the  i  species. 

To  solve  the  differential  equations  a  fourth  order  predictor-corrector  method, 
sometimes  called  the  Modified  Adams-Bashforth  Method,  is  used.  This  method  is 
generally  available  for  computer  use  in  an  Avco  library  subroutine  which  was 
partially  reprogrammed  for  purposes  of  the  problem.  The  technique  is  particu¬ 
larly  attractive  because  of  its  computational  efficiency  and  the  relative  ease 
it  affords  in  maintaining  accuracy  requirements.  The  one  major  disadvantage  is 
that  it  is  not  self-starting. 

Before  the  predictor-corrector  algorithm  can  even  be  applied,  values  of  all  the 
concentrations  Ci(t),  »*],•••,  Ne  and  of  their  first  derivatives  jL  C;(t),  c 

dt 

must  be  known  at  each  of  four  points  in  time.  Information  at  the  first  point 
to=0  can  be  obtained  directly  from  the  initial  conditions.  Values  at  three 
additional  equally-spaced  points  tj,  t2,  and  tj  are  obtained  by  a  fourth  order 
Runge-Kutta  integration  process,  the  interval  A t  between  the  points  being  con¬ 
trolled  ly  program  input.  Runge-Kutta  is  particularly  appropriate  as  a  starting 
method  since  results  at  each  step  depend  only  on  conditions  at  the  last  previous 
step.  In  addition,  a  high  degree  of  accuracy  can  be  established  with  a  suffi¬ 
ciently  small  interval.  (However,  Runge-Kutta  would  present  serious  disadvantages 
if  used  over  the  entire  time  period.  A  relatively  large  number  of  derivative 
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calculations  are  required,  and  no  Immediate  estimate  of  the  truncation  error  Is 
available.)  In  the  fourth  order  Runge-Kutta  method  the  derivatives  are  evaluated 
at  four  points  In  the  Interval  and  weighted  to  give  agreement  with  a  Taylor  Series 
expansion  through  the  fourth  order  term.  For  a  comprehensive  discussion  of 
Runge-Kutta  integration  see  chapter  9  in  Mathematical  Methods  for  Digital  Com¬ 
puters  ,  edited  by  Ralston  and  *ilf(6). 

Given  the  starting  values  obtained  from  Runge-Kutta,  the  predictor-corrector 
algorithm  operates  as  follows .  For  each  concentration  C; ,  i  «  1 , . . . ,  Ng3  a  third 
degree  polynomial  ?/*/  approximating d/dtCi is  fitted  through  the  four  available 
derivative  values.  This  polynomial  is  then  Integrated  to  give  an  exprapolated 
or  predicted  concentration 


at  tjm  +  A  t  From  the  set  Ci  (ttfj,,  i  ■  2,  ...»  Ns  approximate  derivatives 

d/dt  Cj  (tj)  are  calculated  directly.  Then  a  new  polynomial  ^  (t)  is  fitted 
through  d/dt  C^tj)  ,  d/dt  Ci(t2),  and  d/dt  6t-  (t3)  ,  and  d/dt  ‘Ci(t4)3  and  interpolated 
or  corrected  concentrations 


are  computed. 

Since  the  truncation  errors  for  the_extrapolation  and  interpolation  formulas  are 
of  the  same  order,  the  difference  Ci(t4 )  -  (t4)  provides  an  estimate  of  the 

actual  truncation  error  incurred  by  the  integration  scheme  in  the  interval  from 
t3  to  t4.  In  fact,  it  can_be  shown  that,  neglecting  higher  order  terms,  the  true 
value  Ci  (t4)  lies  between  C;  (t4)  and  (t4)  .  Again,  see  reference  (6),  chapter 
8  for  a  general  discussion  of  predictor-corrector  techniques. 

In  the  algorithm  the  quantity  JC,*  (t4)  -  Q  (t4)\  is  tested  to  ensure  that  the 
integration  is  sufficiently  accurate.  A  relative  tolerance  specified  through 

", 

program  input,  multiplied  by  ^  Q (0)9  is  taken  as  an  absolute  upper  limit  for 
_  *'-2 

1  Cj  (t 4 )  -  Ci  (t4)\.  If  for  a^y  i  this  limit  is  exceeded,  the  time  interval  At  is 
halved  and  the  integration  retried.  The  process  of  testing,  cutting  the  interval, 
and  reapplying  the  predictor-corrector  step  is  repeated  until  |  C;  (t4)  -  cj  (t4)  is 
sufficiently  small.  Then  the  predicted  values  (t4),  iml,...,N8  are  taken  as  the 
result  of  integration,  and  d/dt  (T;  (t4),  i-i,...,  Na  as  the  set  of  corresponding 
derivatives.  (The  predicted  rather  than  the  corrected  results  are  accepted  to 

?void  computing  the  additional  set  of  derivatives.)  If,  for  each  i,  \Ci(t4)  - 
(t ^) \  is  at  least  one  hundred  times  smaller  than  the  limiting  criterion,  the 
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interval  4 i  will  be  doubled  before  the  nest  integration  atep  ia  begun.  The  whole 
proceae  la  then  repeated  over  aucceaalve  intervale,  at  each  atep  the  lateet  A« 
and  condltione  at  the  laat  four  polnta  entering  the  computation,  until  the  end  of 
the  time  period  is  reached. 

P>e  program  waa  coded  in  FORTRAN  IV  for  uae  on  the  IBM  360  aerlea  of  coaputera 
Execution  times  for  typical  problene  on  the  360  model  75  computer  range  betveexi 
one  and  two  minutes.  8  n 


* 
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TABLE  I 


FLUOROCARBON  REACTIONS  AND  RATE  CONSTANTS0 


R1  C?3I  +  M  -  CF3  +  I  +  M 

-  2.27  x  1030  T-4-0  exp(-57385/RT) 

R2  C2F6  +  M  -  CF3  +  CF3  +  M 

k2  -  8.40  x  1020  T0*5  exp(-76500/RT) 

R2*  CF3  +  CF3  +  M  -*  C2F6  +  M 

k{  -  7.14  x  1017  T0*5 

R3  CF3  +  M  -  CF2  +  F  +  M 

k3  -  T.57  x  1049  T  "9-04  exp (-92254/RT) 

R3*  CF2  +  F  +  M  -*  CF3  +  M 

k3'  *  1.49  x  1046  T  “9‘04  exp (-2287/RT) 

R4  CF^  +  M  -  CF3  +  F  +  M 

k4  -  6.15  x  10 34  T  ”4-64  exp(-122421/RT) 

R4*  CF,  +  F  +  M  -  CF.  +  M 

3  4 

t4'  »  9.79  x  1031  T-4-64  exp(-2849/RT) 

R5  CF2  +  M  -  CF  +  F  +  M 

k5  -  4.20  x  1026  T2'85  exp (-106000/RX) 

R5'  CF  +  F  +  M  -  CF2  +  M 

lj'  -  6.57  x  1026  T-2-85 

M  |  s  collision  partner  (Argon) 

units  in  calories,  cubic  centimeter,  degree  Kelvin,  mole,  second. 
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TABLE  2 


EXAMPLES  OF  CF3!  AND  C2F6  SHOCK  TUBE  EXPERIMENTS 


Experiment 

No. 

Shock 
Velocity 
( jm/n  sec) 

Chock 

Temp. 

<°K) 

P21 

Total 

Gas  Conc'n 
(mole/cm^)  10® 

Apparent 

Rate  Constant 
(k^cm^/mole  sec)  . 

1:100  CF 

^l  -  Argon  Data 

1 

1.322 

1804 

3.66 

7.94 

2.28  x  108 

2 

1.382 

1928 

3.75 

8.16 

9.08  i  108 

3 

1.441 

2060 

3.82 

8.23 

4.95  x  109 

4 

1.478 

2143 

3.87 

8.45 

5.08  x  109 

5 

1.510 

2214 

3.92 

8.56 

7.73  x  109 

1:100  C2 

F^  -  Argor.  Data 

6 

1.347 

1735 

4.00 

10.8 

9.09  x  107 

7 

1.395 

1850 

4.02 

8.72 

5.69  x  108 

8 

1.440 

1960 

4.03 

8.76 

1,20  x  109 

9 

1.456 

2000 

4.04 

8.76 

1.  38  x  109 

10 

1.533 

2205 

4.06 

8.82 

4.88  x  109 

Conditions  for  complete  dissociation  of  CF3I  and  at  shock  front. 


TABLE  3 


EXAMPLES  OF  CF4  SHOCK  TUBE  EXPERIMENTS 


Experiment 

No. 

Shock 

Velocity 

(mm/^sec) 

aShock 

Temp. 

(°K) 

P21 

Gas  Conc'n 
(raole/cnr)10^ 

Apparent 

Rate  Constant 
(k4,cnr/mole  sec) 

1:100  CF^  -  Argon  Data 

11 

1.559 

2410 

3.75 

1.82 

7.96  x  107 

12 

1.623 

2590 

3.78 

1.63 

3.32  x  108 

13 

1.643 

2640 

3.795 

1.42 

5.61  x  ID8 

14 

1.670 

2720 

3.808 

0.518 

1.04  x  109 

15 

1.715 

2850 

3.83 

0.422 

1.89  x  109 

a16 

1.795 

2970 

4.05 

1.09 

4.70  x  109 

2:100 

CF^  *•  Argon  Data 

17 

1.5C7 

2260 

3.81 

1.84 

2.91  x  107 

18 

1.544 

2350 

3.845 

1.76 

6.12  x  107 

19 

1.591 

2460 

3.87 

1.67 

1.34  x  108 

20 

1.623 

2550 

3.89 

1.57 

3.63  x  108 

21 

1.66 

2650 

3.91 

1.39 

8.20  x  108 

CF^  Data 

22 

2.015 

2145 

19.1 

0.515 

5.10  x  106 

23 

2.075 

2248 

19.4 

0,521 

2.42  x  10' 

ste  constant  evaluated  at  AO  percent  decomposition. 
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Figure  1  Difluoromethyfene  ultraviolet  light  absorption 
behind  shock  waves  into 

(a)  1:100  CF3l-Ar,  shock  veloc.  =  1.377  mm/^sec,  Temp  -  1984°  K 

(b)  1:100  C2F£)-Ar/  shock  veloc.  ~  1.456  mm/fisec,  Temp  -  2107°  K 

(c)  1:100  CF4~Ar#  shock  veloc.  -  1.782  mm/^scc.  Temp  =  3050°  K 

at  4  cm  Hg  initial  pressure.  Oscillogram  writing  =  20  // sec/cm 


ELECTRON  CONCENTRATION,  mole  cm' 


5.0  10.  15.  20. 


LABORATORY  TIME ,  MICROSECOND 


Figure  2  Numerical  Check  of  Nonequilibrium  Shock-Tube 

Program  with  Airchemistry  Stream-Tube  Calculation. 
O  ,  ["Rate  of  Ionization  Behind  Shock  Waves  in  Air. 
I.  Experimental  Results,"  Avco/Everett  Research 
Rept.  105,  (September,  I960)]  . 
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SPECIES  CONCENTRATION,  mole  cm 


LABORATORY  TIME,  microseconds 

27-3844 


Figure  3  Chemical  kinetics  shock-tube  calculation  for 

decomposition  of  1:100  CF3l-Argon  gas  mixture 
(Experiment  No.  2).  Q,  spectroscopic 
measurements  of  CF2  concentration.  Error  flag 
reflects  uncertainty  in  absorption  coefficient  used  in 
data  reduction.  Argon  curve  is  (100)  I. 
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SPECIES  CONCENTRATION,  mole  cm’ 


2 


27-3645 


LABORATORY  TIME,  mlc.oseconds 


Figure  4  Chemical  kinetics  shock-tube  calculation  for 

decomposition  of  1:100  C2F5  -  Argon  gas  mixture 
(Experiment  No.  7).  0,  spectroscopic  measure¬ 
ments  of  CF2  concentration. 
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k3  (expt'l),  cm3  mole**  sec* 


SPECIES  CONCENTRATION,  mole  cm' 


27-3649 


l/T  x  I0\  K’1 


Figure  7  Experimental  rate  constant  for  CF4  dissociation 
Solid  iine  is  the  function 

6.15  x  1(P4  t-4.64  exn(-122421/RT) 

cm3  moie'l  sec'l  from  least  squares  fit  of  data 


SPECIES  CONCENTRATION,  mole  cm 


27-3646 

Figure  8  Chemical  kinetics  shock-tube  calculation 
for  decomposition  of  pure  CF4 
(Experiment  No.  22).  Q,  spectroscopic 
measurements  of  CF2  concentration  in 
reaction  mixture, 
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